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Abstract. Elliptic 6j'-symbols first appeared in connection with solvable models of statisti- 
cal mechanics. They include many interesting limit cases, such as quantum 6j'-symbols (or 
g-Racah polynomials) and Wilson's biorthogonal 10W9 functions. We give an elementary con- 
struction of elliptic 6j'-symbols, which immediately implies several of their main properties. 
As a consequence, we obtain a new algebraic interpretation of elliptic 6j-symbols in terms of 
Sklyanin algebra representations. 



1. Introduction 

The classical 6j-symbols were introduced by Racah and Wigner in the early 
1940's |Rac[ IWij . Though they appeared in the context of quantum mechan- 
ics, they are natural objects in the representation theory of SL(2) that can be 
introduced from purely mathematical considerations. Wilson |Wlj realized that 
6j-symbols are orthogonal polynomials, and that they generalize many classical 
systems such as Krawtchouk and Jacobi polynomials. This led Askey and Wilson 
to introduce the more general g-Racah polynomials [AWlJ. 

The g-Racah polynomials belong to the class of basic (or q-) hypergeometric 
series [HTU] . Since the 1980's, there has been a considerable increase of interest 
in this classical subject. One reason for this is relations to solvable models in 
statistical mechanics, and to the related algebraic structures known as quantum 
groups. 

Kirillov and Reshetikhin [KRJ found that g-Racah polynomials appear as 6j- 
symbols of the SL(2) quantum group, or quantum 6j -symbols. We mention that 
in the introduction to the standard reference |CPj . three major applications of 
quantum groups to other fields of mathematics are highlighted. For at least two 
of these, namely, invariants of links and three-manifolds |Tuj , and the relation to 
affine Lie algebras and conformal field theory [EFKJ, quantum 6j-symbols play a 
decisive role. 

The g-Racah polynomials form, together with the closely related Askey- Wilson 
polynomials, the top level of the Askey Scheme of (q-) hypergeometric orthogonal 
polynomials |KSj . One reason for viewing this scheme as complete is Leonard's 
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theorem [Lj, saying that any finite system of orthogonal polynomials with poly- 
nomial duals is a special or degenerate case of the g-Racah polynomials. However, 
if one is willing to pass from orthogonal polynomials to biorthogonal rational func- 
tions, natural extensions of the Askey Scheme do exist. 

One such extension was found by Wilson |W2j , who constructed a system of bi- 
orthogonal rational functions given by io0g (or, more precisely, wWq) basic hyper- 
geometric series. These form a generalization of g-Racah polynomials that seems 
very natural from the viewpoint of special functions; see also RS . 

Another indication that natural generalizations of quantum 6j-symbols exist 
came from statistical mechanics. The solvable models that lead to standard quan- 
tum groups and quantum 6j-symbols appear there as degenerate cases. Typically, 
the most general case of the models involve elliptic functions. In the 1980's, 
Date et al. |DK[ lD] developped a fusion procedure for constructing generalized 6j- 
symbols from i?-matrices of face models. When applied to Baxter's eight-vertex 
SOS model |ABF1 iBa] . this leads to elliptic 6 j -symbols, which include quantum 
6j-symbols as a degenerate case. However, no identification of these objects with 
biorthogonal rational functions was obtained, nor was their nature as generalized 
hypergeometric sums emphasized. 

In the latter direction, Frenkel and Turaev jFTlj found that the trigonometric 
limit case of elliptic 6j-symbols can be written as loWg-series. A further limit 
transition gives rational 6j-symbols. Moreover, in |FT2j it was found that general 
elliptic 6j-symbols may be expressed as elliptic, or modular, hypergeometric series, 
a completely new class of special functions. In spite of their intriguing properties, 
including close relations to elliptic functions and modular forms, such series were 
never considered in "classical" mathematics, but needed physics for their discovery. 
We refer to [GR2J for an introduction to the subject, with further references. 

Frenkel and Turaev seem not to have been aware of the work of Wilson. Spiri- 
donov and Zhedanov SZ Q ISZ2j gave an independent approach to elliptic 6j- 
symbols, showing in particular that they are biorthogonal rational functions, and 
that they coincide with Wilson's functions in the trigonometric limit. More pre- 
cisely, trigonometric 6j-symbols correspond to certain discrete restrictions on the 
parameters of Wilson's functions. Similarly, elliptic 6j-symbols correspond to 
special parameter choices for Spiridonov's and Zhedanov's biorthogonal rational 
functions. We will be concerned with the larger parameter range, although, for 
simplicity, we will use the term "6j-symbol" also in that setting. 

To summarize, we have a scheme (in the sense of Askey) consisting of classical, 
quantum, rational, trigonometric and elliptic 6j-symbols, see Figure [0 Arrows 
indicate limit transitions. We also give the hypergeometric type of the systems. 
(We use Spiridonov's |Spl| more logical notation 12 Vn, see f|5.1|) below, rather 
than as in |FT2j . for the series underlying elliptic 6j-symbols.) 

Note that the discrete part of the Askey Scheme lies below the classical and 
quantum 6j-symbols in Figure We remark that, once we decide to include 
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Figure 1 . The hierarchy of 6j-symbols 

biorthogonal rational functions, many further limit cases exist (including biortho- 
gonal polynomials and orthogonal rational functions). It seems desirable to classify 
all limit cases, along with their continuous relatives. Many known systems (see 
jXTl IW1 IHM1 EMU EM2 IE3 El EU E21 E3 HSj for some candidates) should fit 
into this larger picture. 

The aim of the present work is to give a self-contained and elementary approach 
to 6j-symbols, which works for all five cases. We will show how to obtain many 
of their properties in an elementary fashion, without using quantum groups or 
techniques from statistical mechanics (although the approach is certainly related 
to both). In the exposition we will focus on trigonometric 6j-symbols. We stress 
that this is not because of any essential difficulties with the elliptic case, but since 
we want to emphasize the elementary nature of our approach as much as possible. 

Our main idea comes from the interpretation of Askey-Wilson and g-Racah 
polynomials given in |Rolj ; see (StJ [Z] for related work. The standard definition 
of 6j-symbols involves three-fold tensor products of representations. This works 
equally well in the classical and quantum case. In [Rolj . we gave an interpre- 
tation of g-Racah polynomials involving a single irreducible representation of the 
SL(2) quantum group. On the level of polynomials, this means that g-Racah poly- 
nomials appear as g-analogues of Krawtchouk polynomials rather than of Racah 
polynomials. Realizing the representation using difference operators on a func- 
tion space (sometimes called the coherent state method |Juj). this yields a kind 
of generating function for g-Racah polynomials, see (|2.16|) below. We may now 
forget about the quantum group and use the generating function to recover the 
main properties of g-Racah polynomials. Our aim is to generalize this approach to 
include all 6j-symbols in Figure ^ keeping the underlying quantum group (known 
as the Sklyanin algebra in the most general case) implicit until the final Section 

The plan of the paper is as follows. Section |21 contains preliminaries; in partic- 
ular we explain in some detail the degenerate cases corresponding to Krawtchouk 
polynomials and g-Racah polynomials (or quantum 6j-symbols). In Section |3] 
we generalize this to trigonometric 6j-symbols, and in Section El we sketch the 
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straight-forward extension to elliptic 6j-symbols. Although the main point of the 
paper is to avoid using quantum groups, we give an algebraic interpretation of 
our construction in Section |3 It turns out that elliptic 6j-symbols appear as 
the transition matrix between the solutions of two different generalized eigenvalue 
problems in a finite-dimensional representation of the Sklyanin algebra. 

Acknowledgements: It is a pleasure to dedicate this paper to Richard Askey, 
who has provided help and support in times when it was needed. It is a natural 
continuation of [Rol, which was presented at Askey's 65th birthday meeting in 
1998, and it is directly inspired by his comments on that talk. Furthermore, 
I thank Eric Rains for crucial correspondence on the Sklyanin algebra and Paul 
Terwilliger for many useful discussions. The figures were made using Paul Taylor's 
package diagrams. Finally, I thank Institut Mittag-Leffler for valuable support, 
giving me the time to type up the manuscript. 

2. Preliminaries 
2.1. Notation. We recall the standard notation for shifted factorials 

(a)k = a(a + 1) • • • (a + k — 1), 
(ai, . . . , a n )fc = (ai)fc • ■ ■ (a n ) k , 

for hypergeometric series 

ai,...,a r _ (01, ...,a r ) k x K 
for g-shifted factorials 

(2.1) (a; q) k = (1 - a)(l - aq) ■ ■ ■ (1 - aq^ 1 ), 

(at,..., a n ; q) k = q) k --- (a n ; q) k , 
for g-binomial coefficients 



r- 1 s 



N 
k 

for basic hypergeometric series 

CL~^ j • • • ■ Cliiy 

}i, ■■■,b s ' 
and for very-well-poised series 

r+ iW r (a; b\, ... , b r -2] q, z) = r+ i 



N 



(q;q)k(q; q)N-k 



k=0 



(q,h,...,b s ;q) k 



a,qa,2, —qa*,bi, . . . ,b r _ 2 
, -as, aqjb\, . . . , aq/b r _ 2 

y> l-aq 2k (a,b 1 ,...,b r __ 2 ;q) 
1 — a (q, aq/bi, . . . 



k=0 



q,z 



1 z k . 



r-2] q)k 
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If one of the numerator parameters equals g~ n , with n a non-negative integer, the 
series reduces to a finite sum. We are particularly interested in the terminating 
balanced iqWq, that is, the case when r = 9, the sum is finite, z = q and a 3 q 2 = 
b\ ■ ■ ■ b 7 . The standard reference for all this is |GRlj . 

To write our results in standard notation, some routine computation involving 
g-shifted factorials is necessary. We will not give any details, but we mention that 
all that one needs is the elementary identities 

(2.2a) (a; q) n = (-l^q^a^q^/a; q) n , 

(2.2b) (a; q) n+k = (a; q) n (aq n ; q) k , 

(2.2c) (a; g)^ fc = (-!)*?( V-AO* Mn 



(? 1_n /o; q)k' 



2.2. An extended example: Krawtchouk polynomials. Our guiding exam- 
ple will be Krawtchouk polynomials, arising as matrix elements of SL(2, C). (In- 
cidentally, they also appear as 6j-symbols, namely, of the oscillator algebra [VK, 
§8.6.6].) For later comparison, we recall some fundamental facts on this topic 
jETllVK] . 

Consider the coefficients K l k = K l k (a, b, c, d; N) in 



A? 



(2.3) (ax + b) k (cx + d) N ~ k = K l x \ 



1=0 



where k G {0,1,..., N} and we assume, with no great loss of generality, that 
ad — be = 1. Note that SL(2) acts on polynomials of degree < iV by 



(2.4) p(x) i-> (cx + d) N p 



ax + b 
cx + d 



and that K l k are the matrix elements of this group action in the standard basis of 
monomials. 

Using the binomial theorem, several different expressions for K\ as hypergeo- 
metric sums may be derived. For instance, 

(ax + b) k (cx + d) N - k = Q x + ^ (cx + d)\ (cx + d) N ~ k 



fk\ (N - j\ ^ k _j cm ^ N _ k _j_ m ^ m+ j 



j=0 m=0 
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Thus, writing m = I — j, we obtain 

min(fc,Z) 

(2.5) 



E 

3=0 

N 



k\ (N - j 



j) V 1 ~ J 



b k c l d N - k - l > 



Ft 



h k-3 c l-3 d N-k-l 
~-k,-l 

-N ; be 



in standard hypergeometric notation. 

Note that the expansion problem inverse to (|2.Hj) . 



N 



(2.6) 



^2 K i ( ax + b)\cx + d) 



N-l 



1=0 



is equivalent to the original problem (replace the matrix 
Thus, K{ is given by a similar formula, namely, 



c d > 



by its inverse). 



Ki 



_^\k+l N -k-liJc J 



b k d 2 F 1 



-k,-l 
-N ; be 



Combining (|2.Hj) and (|2.(jjl . we obtain the orthogonality relation 



TV 



(2.7) 



z=o 

E (7) (^)(- i ) i+m « iY_m "^ fc+i c i+m rf / 



X 2^1 



2 Fi 



-AT ; 6 C 

Now let us introduce the standard notation 

— n, —x 1 



m, —I 1 
-iV ; ~Vc 



K n (x; P} N) = 2 F 1 



-N 



P. 



This is a polynomial in x of degree n, known as the Krawtchouk polynomial [Z 
Writing be = —p, ad = 1 — p and t = cx/d, (12. Hj) takes the form 



(2.8) 



p 



N / Art 



which is a well-known generating function for Krawtchouk polynomials. Our ap- 
proach to the 6j-symbols in Figure ^ will be based on generalizing this identity. 
In terms of Krawtchouk polynomials, (|2.7jl takes the form 



( N )p x (l-p) N ~ x K k {x;p,N)K m {x;p,N) = 6 km 



[i-pY 
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For < p < 1, this is an orthogonality relation for a positive measure, namely, the 
binomial distribution on a finite arithmetic progression. That we get a genuine 
orthogonality stems from the fact that the underlying representation is unitarizable 
for the group SU(2). 

Several other interesting properties of Krawtchouk polynomials are immediately 
obtained from (j2.3|) . For instance, one may consider three bases tk, fk, gk, each 
being of the form (ax + b) k (cx + d) , with different a, b, c, d. The transition 
coefficients in 



are then all given by Krawtchouk polynomials, with different parameter p. Clearly, 
they are related by matrix multiplication: 



A' 



(2.9) 



Knm — ^ ] Knk-Kkm- 



k=0 



In the case = gk, one gets back the orthogonality (|2.7|) . 

From the viewpoint of group theory, (J2.9|) corresponds to representing the group 
law (i.e. multiplication of 2 x 2 matrices) in an (iV + l)-dimensional representation. 
This should be quite familiar when N = 1 and we restrict to SO (2), the rotations 
of the plane, obtaining in this way the addition formulas for sine and cosine. Thus, 
(J2.9|) appears as a natural extension of these addition formulas. 

From the hypergeometric viewpoint, (|2.9|) is an instance of Meixner's formula 
M 



E 

k=0 



k\ 



—k, a 



x 



(1 - zY 



-k,b 



c 

+b-c 



1 — z + xz) a (l — Z + xz) 



a, b 



xyz 



[1 — z + xz){l — z + yz) 

-m, c = —N, with m, n, 



More precisely, it is the special case when a = —n, b 
N integers such that < m, n < N . 

Another consequence of ()2.3|) is obtained by exploiting the multiplicative struc- 
ture of the basis vectors. Namely, expanding both sides of 

(2.10) [ax + b) k+j (cx + d) M+N ' k ' ] = (ax + b) k {cx + d) M ~ k {ax + b) ] {cx + d) N ^ 

into monomials gives 



' n x n 



l m n 

or 

(2.11) K l k+j {a, b,c,d;M + N)= K k^^ *>, c, d; M)K?(a, b, c, d; N). 



m-\-n=l 
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In hypergeometric notation, this is 

-I, -k-j. 



(2.12) 



M + N 



1F1 



-M-N 

= E 

m+n=l 
0<m<M 
0<n<N 



M 

m 



— m, —k 
-M 



-n, -J 
-N 



The group-theoretic interpretation of (|2.11|) is the following. Let Vn denote the 
(N + 1) -dimensional irreducible representation of SL(2), realized on the space of 
polynomials as above. Then multiplication of polynomials defines a map Vm ® 
Vn — > Vm+n- The relation ()2.10|) . and thus ()2.11|) . expresses the fact that this 
map is intertwining, that is, commutes with the group action. This immediately 
suggests a non-trivial generalization. Namely, one has the equivalence of represen- 
tations Vm <£> Vn — ®^^ M Vm+n -2s, and one may do the same thing for the 
intertwiners Vm®Vn — ► Vm+n -2s- The corresponding generalization of (|2.12j) has 
additional factors of type 3_F 2 appearing on both sides. From the group-theoretic 
viewpoint, these are Clebsch-Gordan coefficients and, from the viewpoint of spe- 
cial functions, Hahn polynomials, see VK, §8.5.3]. 

Next we point out that (|2.11|) may be iterated to 

(2.13) K l ki+ ... +kn (a, b,c,d;M 1 + --- + M n ) 

= K 1 ( a ' b,c,d;M 1 )--- K™ : (a, b, c, d; M n ), 

m-iH Ym n =l 

where < h^rrii < Mj. This is especially interesting when Mj = 1 for all i. 
Writing the result in hypergeometric form, we get in that case 



iFt 



I . k ]_ • ■ • fa 

—n 



t 



vn\-\ \-m n =l i=l 

0<mi<l 



-1 ^ 



Note that the range of summation may be identified with the /-element subsets L 
of N — {1, . . . , n} (interpreting mj = 1 as i G L). Similarly, (ki, . . . , k n ) labels a 
subset K of N with J2i h elements. Since 



2*1 



-mj, ki 



t 




TYtj^ — kj^ — 1 . 

otherwise, 



the term in the sum is (1 — £)' inis 'L Replacing t with 1 — t, we now have 



(2.14) 



\N\ 
I 



2^1 



-l,-\K\ 
-\N\ 



l-t 



\Lf\K\ 



K C N. 



E " 

LCN, \L\=l 

This (not very deep) identity gives a combinatorial interpretation for Krawtchouk 
polynomials as a generating function for the statistics \L Pi K\ on subsets L of 
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fixed cardinality. We shall see that the appearance of 6j-symbols in statistical 
mechanics is via a generalization of this identity. 

2.3. g-Racah polynomials. In |Rolj . we considered a g-analogue of the above 
set-up, leading to general g-Racah polynomials. The group SL(2) was replaced by 
a quantum group, and the basis vectors x k and (ax + b) k (cx + d) N ~ k by appropriate 
g-shifted products such as 

k N-k 

Y^axq 1 ' + b) J^J (cxqi + d). 

3=0 j=0 

Such bases were interpreted as eigenvectors of Koornwinder's twisted primitive 
elements |K2j . and also as the image of standard basis vectors x k under Babe- 
Ion's vertex-IRF transformations [B] (called generalized group elements in |Rol| ). 
Actually, we focused on the case of infinite-dimensional representations, and only 
mentioned the case of present interest somewhat parenthetically |Rol[ Section 6] . 

To be more precise, the expansion problem that yields quantum 6j-symbols 
(g-Racah polynomials) is 

k N-k N I N-l 

(2.15) \\(axq •' + b) Yl {cxq- 1 +d) = Y,C l k Y[{axq j + (3) J] W + <*)• 

3=0 j=0 1=0 j=0 j=0 

For generic parameter values, the polynomials on the right form a basis for the 
space of polynomials of degree < N, so that the coefficients exist uniquely. For 
the rest of this section we assume that we are in such a generic situation. 
Note that when q = 1, (j2.15J) reduces to 

N 

(ax + b) k (cx + d) N ~ k = Y,°k ( ax + P)\l x + 

1=0 

which is further reduced to (J2.3|) by a change of variables. The expansion (J2.15|) 
is more rigid. After multiplying with a trivial factor and changing parameters, we 
may restrict to the case 

N 

(2.16) (ax; g _1 ) fc (fa;; q^) N -k = ^ C l k (a, b, c, d; N; q) (cx; q)i(dx; q) N -i. 

1=0 

We could dilate x to get rid of one more parameter, but the remaining 7 para- 
meters, counting q, enter in a non-trivial fashion. Indeed, we have 

d(a, b, c, d; N; q) = q^ ^ ^ 



(2.17) 

X 403 



I J q (q l - N c/d; q)i(q- l d/c; q) N . l (q^ N b/ c; q) k 

q~ k , q~ l , q k ~ N b/a, q l ~ N c/d 

q~ N ,c/a,q l - N b/d ;q ' q 

In Section El we will derive a more general identity in an elementary way. 
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Similarly as for (|2.Hjl . we may invert (|2.1fi|) to get the orthogonality relation 

N 

(2.18) 5 km = C l k (a, b, c, d; N; q) Cp(c, d, a, b; N; q' 1 ). 

1=0 

One may verify that ()2.18|) gives the orthogonality of g-Racah polynomials. (If we 
want a positive measure, some conditions on the parameters must be imposed.) 

Note that (|2.1fi|) generalizes the generating function (|2.8j) to the level of g-Racah 
polynomials. This identity was obtained, in a related but not identical context, 
by Koelink and Van der Jeugt [IKV[ Remark 4.11(iii)]. 

The mixture of base q and g _1 in (|2.16|) is crucial. Admittedly, the expansion 
problem 

N 

(2.19) (ax; q) k (bx; q) N - k = Y D k( a > & ' c ' d > N > ?) ( cx > 9)i(dx; Q)n-i 

1=0 

is immediately reduced to ()2.16|) by a change of variables; explicitly, one has 

D l k (a, 6, c, d- N; q) = C l k {aq k -\ bq N - k ~\ c, d; N; q). 
However, the relation 

k 

hm = D fe( a ' & ' c ' d 'i N 'i 9) A m (c, d, a, 6; A^; q) 

1=0 

is not equivalent to f)2. 18j) . It gives a system of biorthogonal rational functions. 
When cd = ah G R it is an orthogonal system, found in an equivalent context by 
Koelink |Ko| Proposition 9.5]; see also |GM| Corollary 4.4] and Remark l4. 21 below. 

Remark 2.1. We conclude the introductory part of the paper with some com- 
ments on the relation to Terwilliger's concept of a Leonard pair; see |T2j and refer- 
ences given there. As was mentioned above, if we let e k = (ax; q~ l ) k {bx; q^^N-k, 
f k = (cx; q) k (dx; q)N-k, then e k and f k appear as eigenbases of certain g-difference 
operators Y\, Y 2 , respectively. It is easy to check that each of these operators acts 
tridiagonally on the eigenbasis of the other, that is, 

Yifk e span{/jfc_i, f k , fk+i}, Y 2 e k e span{e fc _i, e k , e k+1 }. 

Except for a non-degeneracy condition, this is the definition of (Yi,!^) being a 
Leonard pair. Then (|2.17J) means that (Yi, Y 2 ) is a Leonard pair of "g-Racah type" , 
which is the most general kind. This gives a simple model for studying Leonard 
pairs. For instance, the "split decompositions" |Tlj are easily understood in this 
model. A typical split basis between e k and f k would be g k = (ax; q" 1 ) k {dx; q) N-k, 
which interpolates between the two other bases in the sense that 

g k E span{e fc , e k+1 , . . . , e N } n span{/ , fi, ■ ■ ■ , fk}- 

The corresponding split decomposition is then simply Q) k=0 Cg k . More generally, 
we may picture the factors (ax;q~ 1 ) k , (bx;q~ l ) k) (cx;q)]y- k , (dx;q)^- k as being 
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attached to the corners of a tetrahedron, with two edges corresponding to the orig- 
inal Leonard pair and four edges corresponding to different split decompositions. 

3. Trigonometric 6j-symbols 

It is not hard to check that both g-Racah polynomials and Koelink's ortho- 
gonal functions are degenerate cases of Wilson's biorthogonal functions. Thus, if 
one wants to obtain general trigonometric 6j-symbols in a similar way, it seems 
necessary to unify the products (ax;q) k and (ax;^" 1 )^. The correct unification 
turns out to be the Askey- Wilson monomials h k (x; a) = h k (x; a; q), which are the 
natural building blocks of Askey- Wilson polynomials AW2J. They are given by 

fc-i 

h k (x; a) = — axqi + a 2 q 2 ^). 

3=0 

To write this in the notation (J2.1)) one must introduce an auxiliary variable £ 
satisfying 

(3-1) e + r 1 = a:; 

then 

(3.2) h k (x; a) = (of, a£ _1 ; q) k . 

(In the context of Askey-Wilson polynomials one usually dilates a; by a factor 2 
and writes x/2 = cosd, £ = e iG '.) We will need the elementary identities 

(3.3) h k (x; a) = q k ^a 2k h k (x; q^/a), 

(3.4) h k+ i(x; a) = h k (x\ a)hi(x; aq k ). 
It is easy to see that 

(3.5) lira h k (x ft; at) = (ax;q) k , \imt 2k h k (x/t; a/t) = a 2fc g fc(fc_1) (x/a; q~ l ) k - 

Thus, we may unify (j2.16|) and (|2.19jh together with several related expansion 
problems (see Remark 14.21 below) . into 

TV 

(3.6) h k (x; a)h N - k (x; b) = } j R l k (a, b, c, d; N; q) hi(x; c)h N ^i(x; d). 

1=0 

We will suppress parameters when convenient, writing 

R[ = R[(a } 6, c, d; N) = R l k (a } 6, c, d; N; q). 

Note that, in contrast to the limit cases considered above, we cannot get rid of 
any parameters by scaling x. We shall see that R l k depends on all 8 parameters 
(counting q) in a non-trivial fashion. 
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Clearly, the coefficients R l k exist uniquely if and only if (h k (x; c)fiN-k{ x 'i d))k=o 
form a basis for the space of polynomials of degree < N. Although it is not quite 
necessary for our purposes (see Remark [3. 2j) . we will first settle this question. 

Lemma 3.1. The polynomials (h k (x; c)hN- k (x; d))k=o f orm a basis for the space 
of polynomials of degree at most N if and only if none of the following conditions 
are satisfied: 

(3.7a) c/de{ q ^ N , q 2 - N ,...,q N - 1 }, 

(3.7b) cde{q l -\q 2 ~ N ,...,l}, 
(3.7c) c = d = 0. 

Proof. If c/d = qi with 1 — N < j < 0, then all the polynomials have the common 
zero x — d + d~ l , so they cannot form a basis. Similarly, if c/d — g J with 
< j < N — 1 then x = c + c _1 is a common zero, and if (j3.7b|) holds then 
both x = c + c _1 and x = d + d' 1 are common zeroes. In the case (|3.7c|) . all the 
polynomials equal 1 and clearly do not form a basis. 

Conversely, assume that none of the conditions ()3.7|) hold. We need to show 
that any linear relation 

N 

(3.8) ^ A fc h k (x; c)h N - k (x; d) = 

fc=0 

is trivial. By symmetry, we may assume c^O. Choosing x = c + c~ l in ()3.8|) gives 

Xo(dc, d/c; q) N = 0. 

Since (dc,d/c;q)N = only if (|3.7a|) or (|3.7b|) holds, we have A = 0. We may 
then divide (|3.8j) with 1 — cx + c 2 , giving 

TV 

^ A fc h k -i(x; cq)h N - k (x; d) = 0. 

k=l 

By iteration (choosing x = cq + (eg) -1 in the next step) or by induction of N, we 
conclude that Aj = for all i, and thus that the polynomials form a basis. □ 

We now turn to the problem of computing the coefficients R l k . Recall that our 
derivation of ()2.5|) consisted in applying the binomial theorem twice. The same 
proof should be applicable to ([3.6)1 . once we have a generalized binomial theorem 
of the form 

N 

(3.9) h N (x; a) = 22 C k ( a > b, c) h k (x; b)h N ^ k (x; c). 

fc=0 



AN ELEMENTARY APPROACH TO 6j"-SYMBOLS 13 

In fact, ()3.9J) is solved by one of the most fundamental results on basic hyper- 
geometric series: Jackson's %W 7 summation |GRlj D]. Since we have promised to 
give a self-contained treatment, we give a straight-forward proof, motivated by 
our present view of ()3.9|) as an extension of the binomial theorem. 

We will follow the standard inductive proof of the binomial theorem based on 
Pascal's triangle. First we write 

1in+i(x; a) = h^(x; a)(l — aq N x + a 2 q 2N ). 

To get a recurrence for Cj^, we must split the factor 1 — axq N + a 2 q 2N into parts 
that attach to the right-hand side of ([3.9)1 . that is, as 

(3.10) 1 - aq N x + a 2 q 2N = A k (l - bq k x + b 2 q 2k ) + B k (l - cq N ~ k x + c 2 q 2{N ' k) ). 
We compute 

(1 -acq 2N - k )(l -aq k /c) 
k ~ (l-bcq N )(l-bq 2k - N /cy 
1 ' } (l-abq N+k ){l-aq N - k /b) 

k ~ (l-bcq N )(l-cq N - 2k /b) ' 

assuming that the denominators are non-zero. For the elliptic extension discussed 
in Section 03 it is important to note that this uses the elementary identity 

(3.12) - (1 - xy)(l - x/y)(l - uv)(l - u/v) 
x 

= (1 — ux){l — u/x)(l — vy)(l — v/y) — (1 — uy)(l — u/y)(l — vx)(l — v/x), 
with 

(u,v,x,y) i-> (cq N - k ,bq k ,aq N ,£). 
Combining (|3.9|) and (|3.1Uj) yields the generalized Pascal triangle 

(3-13) CS^ = B h CS + A h - 1 C2L l , 

with boundary conditions 

n° — 1 n N — n N — n 
Iterating (|3.13j) . one quickly guesses that 



(3.14) Cf = q k ^ 



k 



(a/ c, q N k ac; q)k{a/b, q k ab; q) N - k 



q (q k N b/c; q) k (q~ k c/b; q) N - k (bc; q) N ' 

To verify the guess, we plug ()3.14|) into (|3.13J) . After cancelling common factors, 
we are left with 

q k - N -\l - q N )(l - q N ab){\ - q N ac)(l - q N+1 - 2k c/b) 

= (1 - q k )(l - q k ^ab)(l - q 2N ~ k+1 ac) (1 - q^cjb) 

- (1 - ^-^(l - q N+k ab)(l - q N - k ac)(l - q N+1 - k c/b), 
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which is another instance of ()3.12j) . this time with 

(u,v,x,y) i — > (q N+ i\fac, q'^^fac, q N ~ k+ ^ ^fac, q k ~^b\/ a/c). 
This shows that, for generic parameters, ()3.9|) holds with the coefficients given by 

Remark 3.2. Note that our proof did not use Lemma 13.11 We see from the 
computation that the expansion exists uniquely unless there are zeroes in the 
denominators of (j3.11|) . which happens precisely if 

b/ce {q 1 ^ ,q 2 ~ N ,...,q N ~ 1 }, bee {1, g" 1 , . . . , g 1 ""} or b = c = 0. 

As expected, this corresponds exactly to the conditions (|3.7j) . 

Remark 3.3. Plugging (j3.14|) into ()3.9|) and rewriting the result in standard 
notation gives 

t XT ( -Nm -n i-N i i hc ut-i \ {cb,c/b,aE > ,aE > - 1 ;q) N 
s W 7 (q b e; q , q foe, a/c, 6f , b£ ;q,q) = —r — . t . ■ 

{ab,a/b,ct,c£ L ;q) N 

This is Jackson's summation. Essentially the same method was used in [Ro2j 
to obtain extensions of Jackson's summation to multiple elliptic hypergeometric 
series related to the root systems A n and D n . 

We may now compute the coefficients R l k in ()3.6|) by applying the "binomial 
theorem" ()3.9j) twice. For instance, using (|3.4jl we may write 

k 

h k (x; a)h N - k (x] b) = ) j C*(a, c, bq N ~ h ) hj(x; c)h N ^j(x; b) 

k N-j 

= Yl C j( a > C ' b( l N ~ k ) C m~ j ( b i C< l j i d ) h i+m{x\ c)h N _^ m (x] d). 

j=0 m=0 

This gives 

min(fc,i) 

4= E C*(a,c,bq N - k )C?_- j (b,cqi,d). 

3=0 

Plugging in the expressions from (|3.14|) and rewriting the result in standard form 
one finds that the sum is a balanced iqWq series. 

Theorem 3.4. For generic values of the parameters, the coefficients R k in ()3.6|) 

exist uniquely and are given by 

R l k {a,b, c,d; N; q) 



q l(l-N) 



'N 
I 



ac, a/ c; q) k (q bd, bj d; q)i(b/c; q) N - k (b/c; q) N -i(bc; q) N ~ k 



(q l N c/d;q)i(q l d/c; q) N _i(cd; q) N (b/c; q) N (bc; q)i 
x 10 W 9 (q- N c/b; q~ k , q~\ q k ~ N a/b, q l ~ N c/d, cd, q^/ab, qc/b; q, q). 
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Remark 3.5. The special case d = of Theorem 13 . 41 was recently obtained by 
Ismail and Stanton jlHl Theorem 3.1] using different methods. 

Remark 3.6. From their definition, it is clear that R l k have the symmetries 

(3.15a) R l k (a, b, c, d; N) = R l N _ k (b, a, c, d; N) = R k ~ l (a, b, d, c; N), 

and from ()3.3|) we have moreover that 

(3.15b) R l k (a, 6, c, d; N) = g~ 2 © a~ 2k 'R l k {q 1 ~ k / 'a, 6, c, d; N). 

Combining these symmetries with Theorem 13.41 gives further expressions for R l k as 



io Wg sums. These are related via Bailey's classical 10W9 transformations |GRlj . 



On the other hand, the explicit expression in Theorem 13.41 implies many symme- 
tries for R l k that are not obvious from the definition. 

4. Elementary properties 
4.1. Biorthogonality. It is clear from ()3.6|) that the coefficients R\ satisfy 

N 

(4.1) 8 nm = R n(a, b, c, d; N; q) ^(c, d, a, b; N; q). 

We will now show that (|4.1|) gives a system of biorthogonal rational functions, 
which is identical to the one obtained by Wilson |W2j . 

To facilitate comparison with Wilson's result, we rewrite (|4.1|) in terms of the 
functions 

qk{N - k) {g - N . q)n{q k-N c/d ^ bc . q ) h ( q -k d / Cj U- q) N . k (cd; q) N 



and 

S m (/x(fc)) 



where 



{cd) n [ N k ] {b/d; q) k {bc, bd; q) N „ n (b/c; q) N - k 

x R k (a, b, c, d; N; q) 

(q~ N , ac, q 1 ~ N /bd, a/ c; q) n 
{q^ N c/b ]q ) n 

X 10 W 9 (q~ N c/b; q~ n , q n ~ N a/b, q~\ q^c/d, cd, q^/ab, cq/b; q, q) 

q m(N- m) (g-JV ac? ad) q m-N a/h . q) m ( q -™b/a; q) N _ m (ab; q) N 
(ab) m [^] (ac, c/a; q) k {ad, d/a; q) N -k 

x -R™(c, d, a, b] N; q) 
(q~ N , ac, q 1 ~ N /bd, d/b; q) m 
(q l ~ N a/d; q) m 

x w Wg(q- N a/d; q~ m , q m ~ N ajb, q~ k , q k ~ N c/d, ab, q x ~ N jcd, aq/d; q, q), 
fjL(k) = q~ k + q k ~ N c/d. 
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Note that R n has the form 

" j ~ l 1 - gV(fe) + q 2t ~ N c/d 



_ q^^d/b + q^+^cd/b 2 ' 

j — o t — o 

with <jj independent of k, and is thus a rational function in fi(k) of degree n/n. 
Similarly, 

m j ~ l 1 - q'lifr) + q 2t ~ N c/d 



S m (v(k)) = $> III 



q t + 1 fx(k)a/c + q 2t+2 - N a 2 /cd' 

j — t — 



In terms of these functions, (|4.1j) takes the form 



iV 



(4.2) ^2w k R n (ii(k))S m (ii(k)) = C n 5 r , 



k=0 

where 



1 - q 2k - N c/d {q~ N c/d, q~ N , ac, q 1 ~ N /bd, b/d, c/a; q) k k 

and 



Wk l-q~ N c/d {q,qc/d,q 1 - N /ad,bc,q 1 - N c/b,q 1 - N a/d;q) k q 



_ (ba,b/a,dc,d/c;q) N 
n (bc,b/c,da,d/a;q)N 

\ — q~ N a/b (q, q~ N , ac, ad, q 1 ~ N /bc, q 1-N /bd, aq/b; q) n _ n 
X 1 - q 2n ~ N a/b (q~ N a/b; q) n 9 ' 

Thus, we have indeed a system of biorthogonal rational functions. 

We now compare this result with the work of Wilson W2.J, who used the notation 

(z + z _1 
;a,b, c,d,e,f;q 

(ab,ac,ad,l/af;q) n _ n 

= — r io Wg(a/e; az, a z, q/be, q/ce, q/de, q ef,q ;q,q), 

(aq/e;q) n 

where 

abode f = q. 

The normalization is chosen so as to make r n symmetric in a, b, c, d. Assuming 
ab = q~ N with N a non-negative integer, Wilson obtained the biorthogonality 
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relation 

N 



i a a\ \ V ( aq k + a l q k , \ 

(4.4) w k r n y ; a, 6, c, d, e, /; qj 



" x r m [ ^ ' a i ^' c -> d-, f -i q ) ( <%,,„■ 



where 

1 — a 2 q 2k (a 2 ,ab,ac,ad,ae,af;q) k k 
1 — a 2 (q,aq/b,aq/c,aq/d } aq/e,aq/f;q) k ^ 



and 



(a 2 q,q/cd,q/ce,q/de;q)N(q,q n /ef,ab,ac,ad,bc,bd,cd;q) n 
{aq/c,aq/d,aq/e,bf;q) N {q/ef;q) 2n Q 



(in |W2j . the factor g _n and the exponent 2 in a 2 q are missing in the expression 
for C n ). Note that the case m = n = of (|4.4|) is the Jackson sum. 

It is now easy to check that (|4.2|) and ()4.4j) are equivalent. The explicit corre- 
spondence of parameters is 

(a, 6, c,d,e,f) 

N / JV / N f i N 1 N i N I 

I— > (g~ T yc/d,q~~5 \/d/c,q^av cd,q ~^ /bV cd,q^b/V cd,q^V cd/a) 

(which is consistent with the relations ab = q~ N , abode f = q) or, conversely, 

(4.5) (a, b, c, d) i-> (\Zc//, q/dy/cf, ay/cf, by/cf). 

Remark 4.1. Continuous biorthogonality measures for the function r n (not as- 
suming ab = q~ N ) were obtained by Rahman |R21 IR4j . see |Sp3| for the elliptic 
case. 

Remark 4.2. Note that, in view of the limit relations ()3.5j) . any one of the sixteen 
expansion problems 

N 

(ax; q^kibx; q^N-k = (esc; q^iidx; q^N-h 

«=o 

with all possible choices of ±, may be obtained as a degenerate case of ([3.6)1 . It 
is easy to see from Theorem 13.41 that the coefficients C l k are always given by 4^3 
or (equivalently, in view of Watson's transformation |GRlj ) sums. Gupta 
and Masson |GM| worked out all such degenerate cases of Wilson's biorthogonal 
rational functions, finding five different systems. The system in |GM| Corollary 
4.2] is related to the expansion 



N 

! ax; q) k {bx; q~ l ) N -h =/ J C l k (cx; q)i(dx; q'^N-i, 

1=0 
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the system in |GM) Corollary 4.3] to 

TV 

(ax; q) k (bx; q'^N-k = ^Cfc (ge; <?)iV-4, 

«=o 

the system in |GM1 Corollary 4.4] is essentially Koelink's functions [Ko , related 
to 

TV 

(ax; q)k(bx; q) N ~ k = y~]C l k (cx; q)i(dx; q) N - h 

1=0 

the system in |GM) Corollary 4.5] to 

TV 

(ax; q) k (bx; q) N - k = C l k (cx; q)i(dx; q~ x ) N - h 

1=0 

and the system in |GM) Corollary 4.6] is the g-Racah polynomials, related to 

TV 

(ax; g~ 1 ) fc (6x; q~ l ) N _ k = C\ (cx; q)i(dx; q) N -i- 

1=0 

(For some of the systems Gupta and Masson gave a more general version, with 
infinite discrete biorthogonality measure.) All other cases may be reduced to one 
of those five. 

4.2. Addition formula. By iterating ([3.6)1 . one immediately generalizes the bi- 
orthogonality relation (j4.1j) to 

TV 

(4.6) iC(a, b, e, /; N; q) = ^ R k n (a, b, c, d; N; q) R™(c, d, e, /; N; q). 

k=0 

This is an extension of the addition formula ([2.9)1 . We do not believe that the 
general case of ()4.6)) can be found in the literature, although it can probably be 
obtained by analytic continuation from the Yang-Baxter equation for trigono- 
metric 6j-symbols [DKj IFTlj . Though in the present approach it seems almost 
trivial, in a more direct approach, such as defining R™ through the explicit ex- 
pression in Theorem 13.41 it might not be easy to guess the existence of such an 
identity, nor to give a proof. 

It may be of interest to rewrite (J4.6)) in Wilson's notation (|4.3)1 . We introduce 
s = e/a, t = b/f as new parameters, and then make the change of variables 
(|4.5jl . The calculations are essentially the same as those in Section 14.11 and we 
are content with stating the end result. 
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Corollary 4.3. For abode f = q, ab = q N and s and t arbitrary, Wilson's func- 
tions (|4.3|) satisfy the addition formula 



(4.7) ^w k r n i ;a,b,c,d,e,f;q\ 



; a, b, cs, dt, f / s,e/t;q 



X KiVcJf, q/d^cj, s^cjf, q/td^cj; q, N) 
Yr n ( Agm+ f ^ m ;A,B,C,D,E,F;q 



where 



l — aq k (a ,ab,acs,ad,ae,af / s;q) k k 

1 — a 2 (q, aq/b,aq/cs, aq/d, aq/e, aqs/ f;q) k ^ ' 

(a 2 q,q/cdst,qt/ces,q/de;q) N 

X = - — — (ab, ad, bd; q) n 

(aq/ cs, aq/d, aq/e, bj / s; q) n 

v (g, q m st/ef, acs, bcs, cdst; q) m 2m _ N 2_ n 2 2 
(qst/ef;q) 2m 



Y 



(a 2 q, q/de, q/cds, q/ces; q) # (ad, bd; q) n 
{aq/cs, aq/d, aq/e, bf / s; q) N [q/ces, dt/e; q) n 
(ab, t, acs, bcs, dt/e; q) m 



(qs/df,qs/ef;q) m 

(A, B, C, D, E, F) = (y^st/ej, ab^eJJTt, c^e~s~/Jt, d^ft/e~s, ^/ejtfs, ^fejsjt). 

The intermediate expression in ()4.7j) makes it clear that the special case s = 
t = 1 gives back (|4.4j) . since then i?™ = S nm . The presence of square roots is due 
to Wilson's choice of parametrization. Writing the identity explicitly in terms of 
loWg-series, all square roots combine or cancel. 

4.3. Convolution formulas. Next we extend ()2.11|) to the present setting, by 
exploiting the multiplicative property ()3.4|) of our basis elements. Because of the 
shifts appearing in that identity there are several different convolution formulas, 
which we write compactly as follows. 

Corollary 4.4. The coefficients R™ satisfy the convolution formulas 
(4.8) R{ +] (a,b,c,d;M + N;q)= £ Rt(aq a \bq^ N ^ ,c,d; M;q) 

m+n=l 

x R](aq {1 ~ a)k , b q (i-P)(M-k) ^ cg rn dq M-m. N . g j 
for all a, 13 G {0, 1}, where < k,m < M , < j,n < N . 
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Proof. Since, for generic parameters, R l k is determined by (fH.fi)) . it suffices to com- 
pute 

M+N 

^ Q hi(x; c)h M +N-i(x; d), 
1=0 

where Ci is the right-hand side of (|4.8j) . Inside the summation sign, we split the 
factors as 

h m+n (x; c)h M +N-m-n(x; d) = h m (x; c)h M ~ m (x; d)h n (x; cq m )h N - n (x; dq M ' m ). 

Performing the summation, using (|3.6|) . gives 

h k (x; aq aj )h M - k (x; bq^ N ^)h,(x; aq^ k )h N .,(x; bq^^). 

For any a, (3 G {0, 1}, these factors combine to 

h k+j (x; a)h M +N-k-j{x\ b), 

which completes the proof. □ 

4.4. Combinatorial formulas. To get analogues of (|2.1H|) and (j2.14j) . we first 
consider all possible extensions of (|H.4jl to a general sum h kl+ ... +kn (x; a). These 
are naturally labelled by permutations a of {1, . . . , n}: 

h kl+ ... +kn (x; a) = h kam (x; a)h k<a) (x; aq k ^) ■ ■ ■ h K{n) (x; a q k -W + - +k 'l»- 1 >). 

Replacing a by a -1 , this may be written 

n 

h kl+ ... +kn (x; a) = \\h ki (x]aq^), 

i=i 

where we introduced the notation 

{j;a(j)<a(i)} 

for a multi-index k. Note that 

(4.9) \k\? = k 1 + k 2 + --- + k i -i. 

Thus, we have an extension of (|4.8jl labelled by two permutations a, r: 

(4.10) R l kl+ ... +kn (a, b,c,d;M 1 + --- + M n ) 

n 

= E II ^C(«? W S bQ lM ~ klJ , cq^, dq^V; Mi), 

m\-\ \-m n =l i=l 

where < ki,rrii < Mi. (We could replace both occurrences of id in (|4.10j) by a 
third permutation A, but the resulting identity is immediately reduced to (|4.10J) 
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by permuting the m;.) In particular, when Mi 
(4.11) R l kl+ ... +kn (a,b,c,d;n) 



M n = 1, one has 



dq 



\l-m\f 



miH Ym n =l i=l 

0<m 4 <l 



Note that on the right-hand side of (j4.11j) . only the elementary coefficients 



given by 



R™ = R k n (a,b,c,d-l) 



/(l-bc)(l-b/c) (l-bd)(l-b/d)\ 



Rl 



R^ R^ 



(1 
(1 



dc)(l-d/c) 
ac)(l — a/c) 



(l-cd)(l-c/d) 
(1 — ad)(l — a/d) 



\{l-dc){l-d/c) {l-cd){l-c/d)J 

appear. We shall see in Section lo"2"l that the equation (|4.11j) is closely related to the 
fusion of i?-matrices developped in |DKt ID] . This explains the relation between 
our construction and the statistical mechanics approach. 

The combinatorics of the sum ()4.11|) deserves a separate study, but we will 
make some further comments here. Note that, in ()4.11|) . a large number of right- 
hand sides give the same left-hand side. If we only strive for a combinatorial 
understanding of the coefficients R l k , it may be enough to choose the right-hand 
side in a particularly simple fashion. For instance, we may take a = r = id, and 
choose kj as 



(4.12) 



(k 



i- 



k r, 



k n—k 

It seems natural to identify the summation-indices m with lattice paths starting 
at (0, 0) and going right at step i if rrii = 1 and up if rrii = 0, thus ending at (/, n—l). 
Suppose that the i:th step in the path starts at (x,y). Then, by ()4.9|) . 



id 



X, 



lid 



y- 



\1 — m\ 

Moreover, if the ki are chosen as in (|4.12|) . then 

1 = x + y, 1 < i < k, 



k 



id 



k. 



0. 



k + 1 < i < n, 



Ki <k, 



k, k + 1 < i < n. 



U-^lj d 

y 

Thus, for instance, any one of the first k steps in the path that goes right con- 
tributes a factor 

(1 - <f +2 %d)(l - q x a/d) 



R\{aq x+y ,b, cq x ,dq y ;l) 



[1 - q x+ ycd)(l - q x ~yc/d) 
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to the sum. There are three other types of steps, giving rise to similar factors. 
After replacing n by N, this yields the following result. 

Corollary 4.5. The coefficient R l k (a,b,c,d; N;q) is given by the combinatorial 
formula 



V 17 {^~q x+2y ad){l-q x a/d) yr (1 - g 2x+ %c)(l - q y a/c) 
^ 11 (1 - q x +ycd)(l - q x -yc/d) ^ (l - q x +vcd)(l - qy~ x d/c) 

paths early right v * ' y * 11 early up v ^ ' K * 11 



>< n 

late right 



(1 — q x+2y ~ k bd){\ — q x ~ k b/d) (I - q 2x+ y- k bc)(l - q y - k b/c) 



[1 - q x+ ycd)(l - q x - y c/d) 



n 



late up 



q x+ ycd)(l - qy~ x d/c) 



where the sum is over all up-right lattice paths from (0, 0) to (I, N — l), the products 
are over steps in these paths, the first k steps being called "early" and the remaining 
N — k steps being called "late". In each factor, (x,y) denotes the starting point of 
the corresponding step. 

There are many limit cases when Corollary 14.51 takes a simpler form. It might 
be interesting to investigate the limit cases corresponding to various polynomials 
in the Askey Scheme. As an example, let us consider the limit 

L = lim \xm(q k d/bs) N ~ l R l k (as, bs, c, d; N; q). 

s— >0 c— >0 



It is easy to see from Theorem 



I 



that 

(q k - l a/b) k 30i 



q~ k ,q~ l ,q k ~ N a/b q l b 

_-jv > 9) — 

q a 



which, by [GR1, Exercise 1.15] equals 

~N~ 



(4.13) 



(q k a/b) k 302 



q k ,q l ,q k b/a ' 
q -N Q ,q,q 



(As an alternative, one may first use the symmetries (|3.15|) to write 

R l k (a, b, c, d- N) = q-^a-^q-^b-^^Rltf-'/a, q 1+k ~ N /b, c, d), 

and then take the termwise limit in the accordingly transformed version of Theo- 
rem EIU thereby obtaining (J4.13)) directly.) The quantity 1)4.13)1 may be identified 
with a g-Krawtchouk or dual g-Krawtchouk polynomial jKSj . On the other hand, 
starting from Corollary 14.51 gives the combinatorial expression 



e n i n 

paths early right early up 



aq 



k+x 



n in<f=Env n i- 

late right late up paths up early up 

Note that JT q x = g" A ", where || A|| is the number of boxes in the Young diagram 



to the upper left of the path. Writing t = aq ' /b, we conclude that 



'N 
I 



\q- l Mt. 



-N 



q,q 



q llMl t y{k) , 

paths 
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where y(k) is the number of early ups, that is, the ^-coordinate of the end-point 
of the k:th step. This is a simple g-analogue of (|2.14j) . Like (|2.14j) . it is not very 
deep, but it gives an idea about what kind of information is contained in (|4.11|) . 
Note also that when k = or t = 1 we recover the well-known fact 

rivi x - || A |, 

1 paths 



I 



5. Elliptic 6j-symbols 

5.1. Definition and elementary properties. In this section we discuss the 
extension of our approach to elliptic 6j-symbols, or, more precisely, to their con- 
tinuation in the parameters studied in |SZlj . Roughly speaking, this corresponds 
to replacing everywhere "1 — x" with the theta function 

oo 

9(x;p) = -p j x)(l -p> +1 /x), \p\ < 1. 

3=0 

Since 9(x; 0) = 1 — x, the case p = will give back Wilson's functions discussed 
above. The main difference is that in the elliptic case there is no Askey-type 
scheme of degenerate cases; all such limits require p = to make sense. 
We recall the notation |GR2| 

fe-i 

(a;q,p) k = Y[9(aq J ;p), 

j=o 

B{x u ...,x n ;p)= 9(xx,p) ■ ■ ■ 9(x n ; p), 
(ai, ...,a n ; q,p) k = (a x ; q,p) k ■ ■ ■ (a n ; q,p) k . 
Elliptic 6j'-symbols may be expressed in terms of the sum [GR2J 

(5.1) i 2 Vii(a; b, c, d, e, /, g, q~ n ; q,p) 



9(aq 2k ) (a,b,c,d,e,f,g,q n ;q,p) k 



_ . — . . , ; — _ — _ .„ 

^ 9(a) (q,aq/b,aq/c,aq/d,aq/e,aq/f,aq/g,aq n+1 ;q,p) k 

subject to the balancing condition a 3 q n+2 = bcdefg. We mention that this function 
is invariant under a natural action of SL(2,Z) on (g,p)-space, cf. |FT21 Sp2 . 
Since 

(5.2) 9{l/x;p) = -9(x;p)/x, 

the symbols (a;q,p) n satisfy elementary identities similar to (J2.2j) . Moreover, 
(j3.12j) has the elliptic analogue (Riemann's addition formula) 

v 

(5.3) — 9(xy, x/y, uv, u/v;p) = 9(ux, u/x, vy,v/y; p) — 9(uy, u/y, vx,v/x; p). 
x 

As an extension of (|H.2jl we introduce the function 

h k (x;a) = h k (x;a;q,p) = (of, c^" 1 ; q,p) k , £ + f -1 = x. 
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For a ^ 0, this is an entire function of x. (If a = 0, it does not make sense unless 
p = 0.) We may then introduce the coefficients R l k = R l k (a, b, c, d; N; q,p) by 

N 

(5.4) h k (x;a)h N _ k (x;b) = R l k (a, b, c, d; N; q,p) hi(x; c)h N _i(x; d). 

1=0 

Since the computation leading to Theorem 13 . 41 only used results that have verbatim 
elliptic extensions, it immediately carries over to the elliptic case. 

Theorem 5.1. For generic values of the parameters, the coefficients R l k in ()5.4j) 

exist uniquely and are given by 

R l k (a, b, c, d- N ; q,p) = 

(?; q,p)i{q; q,p)N-i 

(ac, o/c; g, p) k {g N - [ bd, bj d; q, p)i(b/c; q, p) N - k (b/c; q, p) N _i(bc; q, p) N - k 
(q l ~ N c/d; q,p)i(q- l d/c; q,p) N ^i(cd; q,p) N (b/c; q,p) N (bc; q,p) t 
x i2V u (q- N c/b; q~\ q~\ q k ~ N a/b, ^c/d, cd, q^/ab, qc/b; q,p). 

Remark 5.2. Like for Theorem 13 A\ the existence and uniqueness falls out of the 
computation, but can also be explained directly. Let / be any function of the form 

N 

(5.5) m = \[B{a j ta i r 1 \p), 

3=1 

and let F be the function 

N 

F(x) = f{e 2wix ) = Y[6{a j e 2wix ,a j e- 2mx ;p). 
i=i 

Then F is an entire function satisfying 

F{x + 1) = F(x), F{x + t) = e- 2mN(2x+T) F(x), F{-x) = F(x), 

where p = e 2mT . In classical terminology, F is an even theta function of order 
2N and zero characteristics. It is known that the space Vjv of such functions has 
dimension N + 1. We will denote by Wfj the space of corresponding functions 
/, that is, of holomorphic functions on C \ {0} such that /(£) = /(£ _1 ) and 
fip€) = (1/p^ 2 ) 7V /(0- ( m jH.alj . these are called BC\ theta functions of degree 
N.) Now we observe that, for ab ^ 0, the functions 

/*(0 = h k (x;a)h N - k {x;b), x = £ + £~ 1 , 

are of the form (|5.5|) . With essentially the same proof as for Lemma f3.1| one may 
check that for 

p m a/b i {q 1 - 1 *, q 2 ~ N , . . . , g^" 1 }, p m ab £ {1, q~\ . . . , q 1 ^}, m G Z, 

(fk)k=o f° rm a basis for Wn- We may then interpret Theorem 15.11 as giving the 
matrix for a change between two such bases. 
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The following Corollary will be used below (in Section 15.21 and in the proof of 
Proposition I6.2J1 . 

Corollary 5.3. If m and n are non-negative integers, then 

h k (x;a)h N - k (x;b) e sp&n k _ m < l < k+n {h l ( y x;aq m )h N _ l ( y x;bq n )}. 

Proof. We need to consider the coefficient R l k (a, b, aq m , bq n ; N; q,p). Isolating the 
factors containing the quotients b/d and a/c in Theorem 15. II gives 

min(fc,£) 

R l k (a,b,c,d;N;q,p) = ^ A i( a / C i Q,P)k-j(b/d; q,p)i-j, 

j=o 

where Xj collects all other factors. If a/c = q~ m this vanishes unless k — j < m, 
and if b/d = q~ n unless I — j < n. Thus, the range of summation is restricted to 
max(/c — m, I — n) < j < (k, I), which is empty unless k — m < I < k + n. This 
completes the proof. □ 

It is clear that the coefficients R k enjoy similar properties as were obtained 
above in the special case p = 0. This applies to the biorthogonality relation (|4.1j) . 
the addition formula (|4.6j) . the convolution formulas in Corollary 14.41 and the 
combinatorial formulas (j4.11|) . In particular, the identity in Corollary 14.41 holds 
after replacing all factors 1 — x with the theta function 9(x;p). 

Remark 5.4. Rains |Ralj has obtained multivariable (Koornwinder-Macdonald- 
type) extensions of the biorthogonal rational functions appearing in Theorem 15.11 
The approach is different from ours, although there are similarities. Note, for in- 
stance, that the one- variable case of the interpolation functions in |Ral[ Definition 
5] are essentially of the form h k (x; a)/h k {x; b). 

5.2. Comparison with statistical mechanics. In this section we compare the 
coefficients R l k with the elliptic 6j-symbols as defined in jDj. We shall see that the 
latter correspond to certain discrete restrictions on the parameters in R l k . 
We will follow the notation of f)KJ, where elliptic 6j-symbols are denoted 

W MN (a,b, c,d\u). 

They depend on four external parameters p, A, £, K and are defined for integers 
a,b,c,d such that 

(5.6) a - b, c-de {— M, 2 — M, . . . , M}, 

(5.7) a-d, b - c G {-N, 2 — N,..., N}. 

As was observed in jFT2j , these symbols may be expressed in terms of the elliptic 
hypergeometric series 12 Vn. Using Theorem 15.11 we may then relate them to the 
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coefficients R k . For instance, we have 
(5.8) WmnU + 2l-N,i + 2k-N, i,j\u) 



(£+j+k+l-N){l-k)+±N(N-u)+±(N-2k)(i-j)(Q_ 



M+M+l-N. 



{q;q,p)n 



R k (a,b, c,d; N; q,p), 



where q = e mX l K and 



(5.9) (a,b,c,d) 



±(u+£+i+l-N) Mu-£.-i+l-N) ±(u+£+j+M+l-N) ±{u-£-j+M+l-N)\ 




or, equivalently, 



(q M , q^ + \ q^ +j , q u ) = (cd/ab, a/6, c/d, abq N ~ l ) 



(£ is a parameter from |DKj that has nothing to do with (j3.1|0 . Note that (|5.7|) 
corresponds to the condition < k, I < N on R l k , while f!5.6|) gives a further 
discrete restriction on the parameters. 

In view of the large symmetry group of the terminating 12^11, there are many 
different ways to identify elliptic 6j-symbols with the coefficients R l k . We have 
chosen the representation (|5.8|) since it explains the relation between fusion of 
i?-matrices and the combinatorial formulas of Section 14.41 Namely, it is straight- 
forward to check that if we let a = r = id in (j4.11|) . specialize the parameters as 
in ()5.9|) and replace i by n + 1 — i in the product, then (j4.11j) reduces to |DKt 
Equation (2.1.21)]. 

It is also interesting to consider the degeneration of the expansion problem 
()5.4|) corresponding to the restriction ()5.6j) . For this we introduce the parameter 
m = (M + j — i)/2. The condition on c — d in (J5.6)) means that m is an integer 
with < m < M. The coefficients (J5.8)) appear in the expansion problem 



1=0 

By Corollary 15.31 R l k vanishes unless k — m<l<M + k — m, which corresponds 
exactly to the condition on a — b in ()5.6|) . Then (|5.10p reduces to 

min(N,M+k— m) 



h k (x;a)h N _ k (x;b) = £ R l k h^x; aq^h^&bq™-" 1 ), < m < M, 



which is thus the expansion problem solved by the elliptic 6j-symbols of [DJ. 

6. Sklyanin algebra and generalized eigenvalue problem 

As was explained in Section 12.31 our approach was motivated by previous work 
on relations between the standard SL(2) quantum group and quantum 6j'-symbols. 
It is natural to ask what "quantum group" is behind the more general case of 



N 



(5.10) 




Z=max(0,fc— m) 
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elliptic 6j-symbols. The answer turns out to be very satisfactory, namely, the 
Sklyanin algebra |ST] . 

We recall that the Sklyanin algebra was obtained from the i?-matrix of the 
eight-vertex model. Baxter found that this model is related to a certain SOS 
(or face) model by a vertex-IRF transformation [BaJ. The original construction 
of elliptic 6j-symbols starts from the i?-matrix of the latter model. Moreover, 
starting from Baxter's SOS model, Felder and Varchenko constructed a dynamical 
quantum group |FV| . which was recently related to elliptic 6j-symbols [KNRJ. We 
summarize these connections in Figure |U1 




SOS model 



D 



8-vertex model 



Dynamical quantum group 



[KNRJ 



Elliptic 6j-symbols 



present 
paper 



Sklyanin algebra 



Figure 2. Connections to quantum groups and solvable models 

It would be interesting to find a direct link between the approach of |KNRj and 
the discussion below. Presumably, this would involve extending Stokman's paper 
|5t] to elliptic quantum groups. In particular, vertex-IRF transformations should 
play an important role. 

To explain the connection with the Sklyanin algebra we introduce the difference 
operators 

A(a,b,c,d)f(0 

r 2 #K, &e, g, p)f(qk) - eo«-\ &T 1 , dg- 1 ; gj/(gjo 

Moreover, N being fixed we write 

A(a, b, c) = A(a, 6, c, q~ N /abc). 

The following observation was communicated to us by Eric Rains, see [Ra2j. 

Proposition 6.1 (Rains, Sklyanin). The operators A(a, b, c) preserve the space 
Wn defined in Remark 15.21 Moreover, they generate a representation of the 
Sklyanin algebra on that space. 
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These representations were found by Sklyanin |S2[ Theorem 4], except that he 
used the equivalent space denoted V N in Remark FT2I and in jSH- Let A i; 

i = 0,1,2,3 be the operators representing Sklyanin's generators Si, pulled over 
from Vjy to W^. Rains observed that every Aj is given by an operator of the 
form A(a, b, c), with specific choices of the parameters and that, conversely, every 
A(a, b, c) may be expressed as a linear combination of the Aj. One may view the 
resulting representation as an elliptic deformation of the group action (j2.4j) . 

Next we consider the action of the operators A on our basis vectors. 

Proposition 6.2. With x = £ + one has 

(6.1) A(a, 6, c)h k (x; q^a)h N _ k (x; q*b) 

q-N 

= —— 9(q k ac, q N k bc, q N ab; p) h k (x; a)h N - k (x; b) 
abc 

and 

(6.2) A(a, b, c)h k (x; A, //) G span^^^^/i^x; q?\, q^fi)}. 

The identity ()6.1)1 is an analogue of the fact that, in the situation of Section 
12.21 any basis ((ax + b) k (cx + d) N ~ k ) k v =Q is the eigenbasis of a Lie algebra element. 
Similarly, ()6.2|) is an analogue of the fact that any other Lie algebra element acts 
tridiagonally on that basis. The parameter shifts are unavoidable and related to 
the fact that elliptic 6j-symbols are biorthogonal rational functions rather than 
orthogonal polynomials, see Remark 16.51 below. 

Remark 6.3. The identities (|6.1|) and (|6.2|) are consistent in view of 

h k (x;a)h N _ k (x;b) G sp&n k _ 1 ^ j ^ k+1 {h j (x; aq)h N -j(x; bq)}, 

which is a special case of Corollary 15.31 

Henceforth we suppress the deformation parameters p, q, thus writing 

9(x) = 9(x;p), (a) k = (a;q,p) k . 

When using notation such as 8(aC, ± ), we will mean #(a£; p)8(a^ 1 ; p) . The follow- 
ing theta function identity will be used in the proof of Proposition 16.21 

Lemma 6.4. If ai ■ ■ ■ a n bi ■ ■ ■ b n+2 = 1, then 

n n+2 
3=1 3=1 



n n+2 

-c^n^^n^r 1 ) 



AN ELEMENTARY APPROACH TO 6j-SYMBOLS 



2!) 



Proof. This is equivalent to the classical identity |TM| p. 34], see also |Ro2j . 

f nyw n h h 

/ . ™ 7T( — / — r — U; ai ■ ■ ■ a n — Oi ■ ■ ■ o n . 

Namely, replace n with n + 2 and bj with bj 1 in that identity, and put a n+ i = £, 
= Moving the last two terms in the sum to the right gives 

a ngw _ riffle nfflcr 1 ^-) 
^K^) egl w " 6(e) nu e W a i) e &- 2 ) ECU ^-v%) ' 

After multiplying with (— l) n £#(£~ 2 ) nj=i a 7 1 $( a j£ ± ) an d using (j5.2|) repeatedly, 
one obtains the desired identity. □ 

Proof of Proposition 16.21 We start with ([6.2)1 . Writing out the left-hand side ex- 
plicitly, collecting common factors and using (|5.2j) repeatedly gives 

(6.3) A(a,b,c,d)((X^ ± ) k (^ ± ) N ^ k ) 

-eOtaC 1 , bC\ c£-\ dC 1 ) (q'hz, qhC^kiq'h^ 



IN-k-l 



-^(ar 1 , &r\ <c\ dc\ q k -hr\ g^r 1 , q N - k ~hr\ g V'r 1 )} • 

Since a&cd = q ,_iV , we may apply the case n = 3 of Lemma 16.41 to the factor 
in brackets. Choose (pi,..., b 5 ) as (a,b,c,d) together with any one of the four 
numbers 

(q k -h,qh- 1 ,q N - k -^,q^- 1 ) 

(let the audience pick it) and choose (01,02,03) as the remaining three of those 
numbers. As a function of £, our expression then takes the form 

(q^^k-^q^^N-k-i 

x {CiB(a 2 (t, o^) + C 2 B(ai^, o^) + C 3 9( ai ^, a 2 f ±) } . 
Depending on the choice of each term is proportional to one of the six functions 

(?~^ ± ) fc (g~5^ ± )j V _ fc , (q^XC^k+^q^C^N-k-i, (?~^ ± ) fc (g^ ± ) JV _ fcj 

(g^A^Or^/i^)^, (g^A£ ± ) fc _i(g - ^/i£ ± )jv-fe+i, (^A^fcCg^^Jv-ifc- 
By Corollary 15.3) these all belong to 

span fc _ 1 < j < fc+1 {(gU^ ± ) j (g5 /i ^ ± ) JV _ j }. 
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This completes the proof of ()fi.2j) . 

If we put A = <pa, /x = q?b in ()6.3j) . the factor can be pulled out 

from the bracket, giving 



A(a,b,c,d) ((qia^kiqzb^N-k 



N-k 



The case n — 1 of Lemma 16 .4| which is equivalent to , now gives (J6.1)) . □ 

Remark 6.5. Proposition 16.21 connects our work with the generalized eigenvalue 
problem (GEVP), which is central to the approach of Spiridonov and Zhedanov 
|SZ1|ISZ2~] . Recall that, roughly speaking, the theory of orthogonal polynomials is 
equivalent to spectral theory of Jacobi operators, that is, to the eigenvalue problem 

Ye k — ^k^k 

for a (possibly infinite) tridiagonal matrix Y. The theory of biorthogonal rational 
functions similarly corresponds to the GEVP 

(6.4) Y x e k = \ k Y 2 e k 

for two tridiagonal matrices Y\, Y 2 . Note that (|6.1j) means that 

e k = h k (x; a)h N _ k (x;b) 

solves the two-parameter family of GEVP:s 

Aie fc = X k A 2 e k , 

where A x = A(q~^a, q~^b, c), A2 = A(q~^a, q~^b, d), with c and d arbitrary. 
Moreover, if we let A3 be any operator of the form A(e,f,g) and we put Y\ = 
A3A1, Y 2 = A 3 A 2 , we have that e k solves ()6.4|) with Yi tridiagonal in the basis 
(efc)^ . Thus, we may view elliptic 6j-symbols as the change of base matrix 
between the solutions of two different GEVP:s, where the involved tridiagonal 
operators are appropriate elements of the Sklyanin algebra, acting in a finite- 
dimensional representation. 
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